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Введение. Транспорт наносов является одним из основных 
процессов, определяющих величины и темпы деформаций 
донных поверхностей водных объектов. Чаще всего прогно- 
стические исследования в этой области строятся на основе 
математических моделей, которые позволяют сократить, а в 
ряде случаев исключить дорогостоящие и опасные в экологи- 
ческом отношении эксперименты. Для прогнозирования из- 
менения рельефа дна в основном используются простран- 
ственно-одномерные модели. Для реальных прибрежных си- 
стем со сложной формой берега вектор потока наносов в об- 
щем случае не ортогонален касательной к береговой линии в 
каждой из ее точек. Также он может не совпадать с вектором 
ветровых напряжений. Поэтому для решения многих практи- 
чески важных задач, связанных с прогнозированием динами- 
ки донной поверхности водоемов, необходимо применение 
пространственно-двумерных моделей транспорта наносов и 
эффективных численных методов их реализации. 

(А. И. Сухинов, 
А. Е. Чистяков, Е. А. Проценко, В. В. Сидорякина) ранее была 


Материалы и методы. — Авторами 
предложена пространственно-двумерная модель транспорта 
наносов, удовлетворяющая основным законам сохранения 
(материального баланса и импульса), которая представляет 
собой квазилинейное уравнение параболического типа. Были 
построены и исследованы линейные разностные схемы и ре- 
шены модельные, а также практические задачи. Однако оста- 
лось в тени теоретическое исследование «близости» решений 
исходной нелинейной начально-краевой и линеаризованной 
непрерывной задач, на основе которой была построена дис- 
кретная модель (разностная схема). Особый интерес пред- 
ставляет исследование корректности линеаризованной задачи 
и определение достаточных условий положительности реше- 
ний, т.к. только положительные решения задачи транспорта 
наносов имеют смысл в рамках рассматриваемых моделей. 
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Результаты исследования. Исследуемая нелинейная двумер- 
ная модель транспорта наносов в прибрежной зоне мелковод- 
ных водоемов учитывает следующие физически значимые 
факторы и параметры: пористость грунта; критическое значе- 
ние касательного напряжения, при котором начинается пере- 
мещение наносов; турбулентный обмен; динамически изменя- 
емая геометрия дна; ветровые течения и трение о дно. Линеа- 
ризация осуществляется на временной сетке — нелинейные 
коэффициенты параболического уравнения берутся с запаз- 
дыванием на один шаг временной сетки. Далее строится це- 
почка взаимосвязанных по начальным условиям — финаль- 
ным решениям цепочки линеаризованных смешанных задач 
Коши на равномерной временной сетке, и таким образом 
осуществляется линеаризация в целом 20 нелинейной моде- 
ли. Ранее авторами были доказаны существование и един- 
ственность решения цепочки линеаризованных задач, получе- 
на априорная оценка близости решения цепочки линеаризо- 
ванных задач к решению исходной нелинейной задачи. В 
данной работе определены условия положительности ее ре- 
шений и их сходимости к решению нелинейной задачи транс- 
порта наносов в норме Гильбертова пространства Г.1 со скоро- 
стью О(т), где т — временной шаг. 

Обсуждение и заключения. Полученные результаты исследо- 
вания пространственно-двумерной нелинейной модели транс- 
порта наносов могут быть использованы при прогнозирова- 
нии нелинейных гидродинамических процессов, повышения 
их точности и надежности в силу наличия новых функцио- 
нальных возможностей учета физически важных факторов, в 
том числе уточнения граничных условий. 


Ключевые слова: пространственно-двумерная модель транс- 
порта наносов, прибрежная зона, нелинейная задача, линеари- 
зованная задача, положительность решения. 
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4. 


Введение. В последние десятилетия количество работ, посвященных рассмотрению нелинейных краевых за- 


дач, заметно возросло. Этому факту способствовало, в частности, появление новых приложений В физике океана и 


атмосферы [1—5]. 


В данной работе рассматривается нелинейная начально-краевая задача для пространственно-двумерной моде- 


ли транспорта наносов мелководных водоемов и проводится анализ некоторых качественных свойств ее решений. 

Исследование указанной задачи осуществляется путем ее линеаризации. С этой целью строится временная 
сетка, на которой значения нелинейных коэффициентов задаются в левом узле соответствующей временной ячейки. 
Далее осуществляется решение цепочки связанных по начальным — конечным данным линейных уравнений парабо- 
лического типа. 

Ранее коллектив авторов проводил исследования линеаризованной начально-краевой задачи, которые вклю- 
чали доказательство существования и единственности ее решения с одновременным определением требований к клас- 
сам гладкости входных данных задачи [6, 7]. Также была получена априорная оценка близости решения цепочки ли- 
неаризованных задач к решению исходной нелинейной задачи в норме Гильбертова пространства Г/1 в зависимости от 
интегральных оценок правой части, граничных условий и нормы начального условия [8]. В предыдущих работах ав- 
торов для линеаризованной задачи транспорта наносов была построена и исследована консервативная устойчивая раз- 
ностная схема, которая численно реализована для модельных и реальных задач прибрежной зоны [9-11]. В данной 
работе определены достаточные условия положительности решений линеаризованной задачи и их сходимости к ре- 
шению исходной нелинейной начально-краевой задачи транспорта наносов в норме Гильбертова пространства Г.1 при 
стремлении к нулю параметра — шага временной сетки, на которой осуществлялась линеаризация. 

1. Постановка нелинейной начально-краевой задачи для двумерной модели транспорта наносов. Следуя 
[4], рассмотрим уравнение транспорта наносов: 
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Здесь Н =Н [% у,1) — глубина водоема; = — пористость донных материалов; 1, — вектор касательного тангенци- 
ального напряжения на дне водоема; т. — критическое значение тангенциального напряжения; т». = аз фу, фи — 
угол естественного откоса грунта в водоеме; А =й{Н ых, У,1) — нелинейный коэффициент, определяемый соотноше- 


нием: 

Аба | те 
А о-в втааН| , 2) 
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гдер1, Ро —— плотности частиц донного материала и ВОДНОЙ среды соответственно; © — ускорение силы тяжести; 6 





К= 


—Щ частота волны; Аи В =—= безразмерные постоянные; а —_ характерные размеры частиц грунта. 


Пусть Рей" — область, где происходит процесс, и 5 — ее граница, которая является кусочно-гладкой ли- 


нией. Областью задания уравнения (1) считаем цилиндр Цт =Рх (0,7 ) высоты Г с основанием ОД. Его граница со- 
стоит из боковой поверхности 5 х [0, Т ] и двух оснований: нижнего Бх{0} и верхнего Бх{Т } . Далее для простоты 


уравнение (1) рассматривается в прямоугольной области 8% У) = {о <х<[,,0<у< в : 


Дополним уравнение (1) начальным условием, предполагая, что функция начальных условий принадлежит 
соответствующему классу гладкости: 











Н (х,у,0) = Но (5), Но(ж5) в С? (Р)оС(Б), втаа(, „Но = С(Ъ), (х,у)ЕБ. (3) 
Сформулируем условия на границе области Р: 

|050 (4) 

Н(Е,,у,#)= Но (у,1), О<у<Г,, (5) 

Н (0, у,#) = Н, (у,#), О<у<Г,, (6) 

Н(х,0,#) = Нз (х), 0<х<Г,, (7) 

Н(х,1//1) = Н. (х,0), бат. (8) 


Дополнительно к граничным условиям (5)—(8) предполагаем выполнение условий их гладкости — существо- 

вание непрерывных производных на границе области р : 
втаа\,)Н © с(ц,) ос (Ц»). (9) 
Считаем, что всегда есть слой ЖИДКОСТИ конечной ТОЛЩИНЫ В рассматриваемой области и для указанного 


временного промежутка не происходит осушения области, т. е. 


Н(х,у,!)> со = соп5Ё > 0, О<х<[,, О<Зу<Е,, 0<1<Т. (10) 





Условие невырожденности оператора задачи имеет вид: 
К > № = сопза > 0, У (х,у)Е РБ, 0<Е5<Т. (11) 


Вектор тангенциального напряжения на дне выражается с использованием единичных ортов системы коорди- 
нат естественным образом: 


=Ин. и Ва в, (ЖИ) (12) 


2. Линеаризация нелинейной начально-краевой задачи транспорта наносов. Для получения линеаризо- 
ванной модели на временном отрезке 0<#<Т построим равномерную сетку ®. с шагом т, т.е. множество точек 


= =, М, МЕТ, 
Осуществим линеаризацию указанной начально-краевой задачи, используя равномерную временную сетку 


т 
0. = {1 =и-т, И=0,1,..., М, М№М-т= т} . Линеаризацию члена Я (коэффициента К) выполним пу- 
т фо 


тем выбора его значений в момент времени 1 =1 


п=0,1,2,..., М-1. 


п=0,1,...,№ и рассмотрения уравнения (1) на временном проме- 


п? 


жутке 1, <Е<&,.1, 
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—1 
При этом предполагается, что известны функция Н (") (ху )=Н й } (х, Ут) и ее частные производные 
по пространственным переменным. 


1 
Если п=| ‚, то в качестве НО (х, У) достаточно взять функцию начального условия, т.е. 


н() (хуи) =Но(х,у). Если же п=2,...,М№ , то функция н“”) (Жув = но (х, у) предполагается известной, 


поскольку предполагается решенной задача (1)—(14) для предыдущего временного промежутка {,„_› <<. 


Введем обозначения: 
В-1 
‚п=1,2,...,М. (13) 


Сре ("-1) 
ть — - эта Н ху, _ 
и фо 8 ( У п 1) 


к) Аба 
(р — Ро ) за) 











Тогда уравнение (1) после линеаризации запишем в виде: 





(1-2). й: > Е до [к га вишн аъ (к 5», <<, ПБ... М (14) 

и дополним его начальными условиями: 
н® (х, у, ю) = Но (х,›),Н() (У) = ИН 1% У, Ва (х, у) Е р, п=2...., М. (15) 
Член вида ак" 5») является при такой линеаризации известной функцией правой части. Граничные 


условия (4)-(8) предполагаются выполненными для всех промежутков времени 1, <#<&,п=1,2,...,М№. 
(1) ы 
Отметим, что коэффициенты А ›, п=1,2,..., М зависят от пространственных переменных х,у и временной 


переменной {„_, п=1,2...., М , определяемой выбором шага т сетки ©. „т. е. к") = ро (У зы п=1,2,... М. 


3. Положительность решений линеаризованной начально-краевой задачи транспорта наносов. Суще- 
ствование и единственность решения линеаризованной задачи (14) с начальными условиями (15) и с условиями на 


границе (4)-(8) обосновано в работах [6-8]. Обозначим решение линеаризованной задачи через Я [% ,1), (% У) ЕД. 


Решение А (х, у. 1), (х, у) Р зависит от шага т сетки ©. ={ = п=0,1,..,№, №-х=Т}. По существу имеем се- 


мейство решений | (") (х, 1] ‚ п=1,2,..., М, зависящее от параметра т. Убедимся в положительности решений ука- 
занной линеаризованной задачи. 
Под положительным решением рассматриваемой начально-краевой задачи понимается функция 


Я“) [УР ЕВ ВМ влавсас” (Ц) ос(Цу ) | ата „В“ Е С(Ц') ‚ удовлетворяющая уравнению 


(14) с начальными и граничными условиями (4)-(8), (15), для которой Я“) (х,у,1) >< >0, су = сопзь, п=1,2,..,М 
при любом и=1,2,...,М№. 


Относительно рассматриваемой задачи будем предполагать выполнение условия 


& к” С | <0 (16) 
для слагаемого левой части уравнения (14), а также положительность начального условия (15): 
Но(х,у),> сд = сопзй > 0, (х»у)ЕР (17) 
и граничных условий (5)-—(8): 
Ни (у, 1), Но (у,1), Нз(у,6), Н(у,1) > в = сот >0. (18) 


Рассмотрим вначале случай и =1, для которого существует единственное решение уравнения (14). Обозна- 


чим данное решение Я () (х, У, 1) . Указанное решение принадлежит классу С ь (Ц й о С (ц. ) ь &та4| „, ЯН ЕС (Ц, ) : 


т (1 
Предположим, что функция Я } (х, у,#) достигает своего минимального значения су или значения меньшего 


его в некоторой внутренней точке параллелепипеда В =В х[0, 1] или в точке его «верхней крышки». Тогда это про- 
тиворечит принципу максимума (минимума): решение начально-краевой задачи не может достигать экстремума, в том 


числе минимума во внутренних точках области В = Бх[0, 1] и на ее «верхней крышке», где во всех точках, кроме 


Сухинов А. И. и др. Достаточные условия сходимости положительных решений 





относящихся к боковой поверхности 5х[0,1 | параллелепипеда В =Бх[0, 1] (рис. 1), ищется решение начально- 


краевой задачи для уравнения (14) при п=1 [12]. 





Рис. 1. Расчетная область 


Е!1о.1. СотриаНопа]! доташ 


т (1 
Следовательно, функция Йй( } (х,у,г) достигает экстремума, в том числе минимума, на нижнем основании 


рх {0} параллелепипеда А = рх [0, 1] или на его боковых границах 5х [0,1 |] „т.е.: 


20 (х, у,г) > шт шт Я (х,у,1), шш я (х,у,1) (19) 
(х.у)еБи=0 (х,у)Е5,0<1<и . 
В силу условий (17), (18), имеем 
тт Я (ху) = >0, шш НО (ху, = >0, 
(х,у)еБи=0 (х,у)е5,0<1<1 
и потому 
я (х, у,Г) > с = сопзй > 0. (20) 
Таким образом, финальное решение Я () (х, У, 1) для 0 <{< 1 будет положительным и большим су. 


2 1 
При п=2 имеем начально-краевую задачу с начальным условием Н ( )(худ)=Н ( } (х, у, ), имеющим 
гладкость такую же, как и начальное условие для уравнения (14) номера и =1. Рассуждая аналогично, решение урав- 


нения (14) 1 <#<ь существует и единственно. Обозначим данное решение Я (2) (х, у,г). Известно, что функция 
я) (х, у,г) принадлежит классу в. (4. )ос(ц, | та4(, Я) в с(ц | 
2 2 х,у) ь 

Очевидно, что решение Я (2) (х, У,1) второй задачи цепочки начально-краевых задач вида (4)—(8), (14), (15) 
достигает своего минимального значения су или значения меньшего со в некоторой внутренней точке параллелепипе- 
да В, = Ррх [4,5] или в точке его «верхней крышки», что противоречит принципу максимума (минимума). Тогда ре- 
шение Я“) (х, у,!) начально-краевой задачи достигает экстремума на нижнем основании указанного параллелепипе- 
да или его боковых границах 5х [а - ь] ‚ и потому выполняется неравенство: 

Н@) (х,у,1) > ©) = соп5 > 0. (21) 
Отсюда следует, что решение Я () (х, у,г) для 0 <Е< и будет положительным и большим величины су. 


Нетрудно видеть, что вообще при любом п=А, К =1,2,..,М№М получаем смешанную задачу для линейного 


уравнения параболического типа, имеющего достаточную гладкость начальных и граничных условий. Для любого 
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номера КА, А=1,2,..,№ функция решения Я (х,у,4), ЗВ, =Ь2...М класса с?(ц,)ос(ц, ) , 


ата „,)Н (ЕС (ц,) начально-краевой задачи (4)—(8), (14), (15) в параллелепипеде Р; = Вы В | достигает экс- 
тремума на нижнем основании указанного параллелепипеда или на его боковых границах, где в силу условия 
(17)-(8) будет выполняться 

Я (х, 7,1) > 6 = сои > 0, К, К=1,2,...,М. (22) 


Учитывая полученные соотношения (20)—(22), приходим к теореме. 
Теорема 1. Пусть даны уравнения (14): 


(п) 
п = к) ве. отааН®) -4ь(«"”) 5}, 1<#<1,П=Т2,..., М, 
01 $1 фо 
в. Неа, Г); Ве. 
В-1 





Аба , т, Тс - тааН (х, У, в) 
(р - ро) ва" Уфо 


н( (х, у, )= Но (ху), н“”) [ЖЕ = н("-)) (хуи), (ху) е р, п=2,...М и с граничными условиями (4)-(8), 


Здесь А) = ь аъ) <0 с начальными условиями 








удовлетворяющими условию (17)—(18). 


Тогда, если  К"О> № > 0, ЕС! (2) | то для любого п=1,2,....№ решение 


Я) (х, ,1), 11 <ЕЁ<&, П=1,2,...М уравнения номера п в цилиндре Цут =Бх(о,Т) является положительным и 
выполняется оценка 


<<, П=Ь2,...,М. (23) 


п 


Я”) (х, у,Е) > су = соп5ё > 0, 1 


4. Сходимость решений линеаризованной задачи к решению нелинейной начально-краевой задачи. 
Убедимся в том, что решение линеаризованной задачи стремится к решению нелинейной задачи в норме Гильбертова 


пространства (р) при ®-—>0, №М-т=Т . С этой целью обозначим решение нелинейной задачи через 
Н(х,у,1), (ху) ЕР, а решение линеаризованной задачи — Н(х,у,#), (х»у)еР. Будем считать, что функция 
Н(х, У,1) ограничена на промежутке 0 <{<Т ‚ а также существуют и ограничены смешанные частные производные 


до третьего порядка включительно по пространственной и временной переменной { от функций Н (х, У,1) ‚ Н (х, У, 1) 


ое} ее, 32 


Заметим, что функция 2 1% 5,1) ограничена на каждом временном промежутке #1 <#<&, п=1,2...., М, по- 


для всех (х,у)еР, 0<{<Т: 








скольку является решением линейной начально-краевой задачи для уравнения параболического типа с «гладкими» 
входными данными. 


Обратимся к нелинейной задаче (1)—(8). При подстановке функции решения Н [1 У,1) в уравнение (1), полу- 


ЧИМ: 





ОН й Пре Е = 
(1-=). — [ет ‘внин шут) (24) 


В-1 





- атааН 


Аба - 
га ры [63] Е _ Тре 


ЕАН. = В 
(р — ро) 4" 51 фо 


Умножим обе части уравнения (24) на функцию Н (% У,1) ? (х, у) ЕР и проинтегрируем по переменной 





1, 0<#<Т ипо переменным (% У) в области О. Получим 
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ры 











(25) 
т т >. 
-| те “ Зе. кеоан фхау @- [ [рн . (+ “ть рана. 
т фо о\ р 
Выполняя ряд и, в равенстве (25), получим: 
5(- =) (и 2 (х,у,Т)-Н? (№ у,0) аа = 
я я (26) 
х Ц 
= | ри. 2 а ‚риа. 
о\р 51т фо о\р 
Далее рассмотрим линеаризованную задачу. При подстановке Я (х, у,г) в уравнение (14) получим: 
ЕТ (п) ы > В 
(1-=).2 ам #9. Че руваН 9 а (Е 3 } НЕ, ПЕЬ,...М, (27) 
0: Ш _. 
Аба м 
где "9 = ы - о ве орааН(" 9 
(р -ро)=4) о 





Умножим обе части уравнения (27) на функцию Я“) = Я") (х,у,!), (х,у)е р и проинтегрируем обе части 
этого уравнения по области Д , а после этого — по переменной & #,_ <Ё<&.. Получим: 


1 


ь >. ЕТ (п) 
| Голод -- о = 
т р 


п- 


_ | [я 4" бани а Г [ря Е" па 
р 31 фо ,:^В 


п 





(28) 


Просуммировав обе части соотношения (28) по п, и =1...., М , находим: 


м6 да ая” 
и-5 ГЕО. - рб - 
р а 


= у | | [я бей _ | [ря Е ла «| 
п=И (р 31 фо р 


т 








(29) 
Преобразуя левую часть равенства (29), запишем: 
5-й? (х,у,Т)-Н? (х, у,0) = 


в | | 4" Вены лы | [ря 45 (Е #") пам 
и Ут фо р 


= В в 


(30) 


Вычитая из равенства (30) выражение (26) и учитывая, что Н (х, у,0) = Я (х,у,0) , получаем: 


5) (и? (х,у,Т)-Н? (х, Тиха = 





Туя (+ не | о а ДП изв }ан (31) 


р (2) 


А ут фо р 


п-1 


| | | п 4 Е" 1). т ‚ тай пы - | [ря аа] 


Выполним преобразования правой части равенства (31). С этой целью добавим, а затем вычтем выражения 


Я) 4 те 
51 





ГИ . — 
} Й( ) . {Е *ТЬ ) под знаком интеграла. Комбинируя слагаемые, получим: 
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5-е) (и (ху,Т)- И? (х,>Т) ри > 





= Ы ея) Чье ‚ таа Н | 
п=1 (В $11 о 
(32) 





+ТЯ(“) «| Тре ‚(к зтаа Вай пана" + 
р зт фо 


Ня“ -Н) «ав (к-ль)акау+ 18" ‚4 (Е) -®).ъ р) «| 


Оценим каждый из интегралов, стоящих в правой части равенства (32) под знаком суммы. Введем обозначе- 





ния: 
п _ |: — АИ”) . Пре Вр 
Ц = Го Н [ет отн о (33) 
А) ут фо 
Б = Ра Чье (к-та Вне а (34) 
В 511 Фо 
та Рае нае) а, (35) 
СР 
П= 1 [ря ‚в (Е -- 5 ура п=1,.... М. (36) 
О 


Полагая п =1 в равенствах (33)-(36), рассмотрим интегралы по первому временному интервалу 1 <Е<4. 
Для начала выполним оценку интеграла /\. Используя разложения функций Н (х,у,г) и я (х, уг) по фор- 
муле Тейлора относительно (1—1) с остаточным членом в форме Лагранжа, с учетом условия (3), можем записать: 


дН (х, у, ) а) (х,у, Е.) (и) 
0+ бе ие 





Н(х,у,^)- Я (ху, 1) - (37) 


ГДе 1 <& <Е, Ц <Ё, <, Ц <ЕЁ СЦ. 
Используя равенство (37), получаем выражение для интеграла / ‚ : 


1 
П= тм 


п 


Й (1) х 
«ее ОН ( о ет] 
р эт фо 


м 0 О 


Л (в тет пан об = 


ур) 511 Фо 
(38) 





С учетом равенства #{—к =т из (38) можно перейти к оценке: 


1.12 1 
И о РМ (39) 


где 





д (х,у,&`)_9Я® (ху, 
м! = тах 3 тах_ е - &) в ы. 2) - У К. отаанН 
<#<й | (х,у)ЕБ Ге [2 зш фо 


Далее получим оценку для интеграла р : Заметим, что справедливо равенство: 


4 . ой. внмн) _ ше в втмй = 
$1 фо $1 фо (40) 


= | ( _ | . и ани) + = ы а. . (зан = стай )} 
$1 фо п фо 


Нетрудно видеть, что интеграл /} с учетом выражения (40) можно представить в виде: 
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П= И я «|5 ((‹ -Е®). втаан + ‚(втааН — втоай ем (41) 


юр о 


Разложим функции К (55 у ) ‚ Н (х, У, 1) по формуле Тейлора и, воспользовавшись равенствами 
К(х, у )= #0) [НЮ я) (х, у, 1% ) ‚ находим: 


=(0) _ О (х, у, Е ) (в) 


КК 42 
9, (42) 








зта4 Н - этаа Я = зта4 Е (ео 4) (2-0 ) (43) 


ГДе 1 <Е 54, Ц < 64 <Е, Ц <Е ЗЕ. 
Используя выражения (42) и (43) для интеграла вида (41), получаем: 


п = и \ УЗКО д [бы отааН + ©). стай В) а (44) 
р 


ут фо 01 


и 


Предполагая ограниченность соответствующих производных, введем постоянную величину: 


- ОК (х, у, : ОН (х, у, 
м} = тах 3 тах_ и вс | ь ый 68) ‚тайн + ©) бб 
ГА 





1 << (х,у)ЕБ Ш фо ОЕ 
Из равенства (44), используя обозначение для величины М}, получим 
1 
1 2 1 
551 М» 1, -1.. (45) 


Для оценки интеграла / - воспользуемся равенством (37), с учетом которого имеем 


1 (1) х И 
н-&-5) м (ель) _оН(н,ве) 





] - > [я деду |, 


где И <5; <&, 11 < <Ь <1<4. 


Вводя обозначение 





1 
М3 = шах $ тах 
10 << (х,у)ЕБ 


ан (х, у, 65) дН (х, у, Ев) а» (къ) 
Ге) 0: ? 


можем записать: 


1 
В т, -Ё, М}. (46) 


Далее перейдем к оценке интеграла 12. С этой целью выполним несложные преобразования /2 с учетом вы- 
ражения вида (42). Получим равенство 
й > ОК (х, у, - 
П= | (1-0). Я. аь. ое аку |. (47) 
и р ОЕ 
В силу ограниченности функции Я (9 (х, ,Г), ® <Е<1 и смешанных производных функции Н (х,у,!) по 


пространственной и временной переменной из (47) получаем неравенство: 


1 
ПЗт М4, (48) 


я [9 


Собирая полученные оценки из (39), (45), (46), (48) и пользуясь неравенством многоугольника для модулей 


где 


<< (х.у)еБ 


м! = тах 3 тах | 








величин, получаем неравенство вида 
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^ 1 * 
Пн" (>. 902 (жу, в) 5.1, -Ё,-М*, (49) 
где 
М! = (м! +М}+М}+М}). 
1-5 
Меняя местами функции Н” (х, ув), Я (ху) и проводя рассуждения, аналогичные предыдущему, 


начиная с равенства (24), можно получить оценку: 


(Я (худ )-Н” [% у, в) <? ЕР -М`. (50) 
р 


Из неравенств (49) и (50) следует неравенство 


ДЕЗ (хуи )-Н? (х, у, )ахау < 2.1, В Му. (51) 
р 





Преобразуем левую часть неравенства (51): 
| | (%у,4)- Н? (ху, в 4 - 
и 7 й (52) 
ы И Я (х,у,н)-Н(х, и) а (%у,6)+Н(ж, и). 
р 
Далее будем предполагать, что выполняются условия теоремы 1 и неравенство вида (23) выполняется для 


функций Н (х,у,!) и НЯ [нету а 


Н(жу >24 >0: (х,у) ЕР, 0<2<Т, (53) 
Я) (х,у,1) >60 >0, (х,у)еБ, 01257. (54) 
С учетом выражений (53), (54) находим 
р (у )-Н? (орлу - 2со р ПА® [5 у, )-Н(х, УВ 2. (55) 
р р 
Из соотношений (51), (52) и (55) следует оценка 
>. 1 * 
ПИ“ (ум) - Н (оз) <—.7 5 Г, "М: В (56) 
Ур) 2 | - 
Требуемая оценка для п =1 получена, т. к. неравенство (56) эквивалентно 
ы 1 
Яхве | <— 42.2 -[.-М*. 57 
ея) ыы ур М (57) 


Фактически сделан первый шаг индукции при и =1 ; к этому относятся выкладки (37)-—(57). 
Пусть и=5, 5=2,.., М. Тогда, рассуждая аналогичным предыдущему образом, получаем оценку: 








— ь: 1 2 2 а 
Я (Т)-Н(Т ЕЕ ЗЕ, С.М (1 + (1-=?)+...+ (1-х ) 58 
| и Л 26 (1-1) 7 ие оо 
Здесь ›‚ -_ Некоторые постоянные функции, определяемые для временного промежутка соответственно. 


Поскольку № = г ‚ то из неравенства (58) получаем: 


[и (т)-Н(Т] < 0-9" -1, с-т-(а-5)+(-2)+..+(1-=“)) (59) 
При >> 0 из (59) находим 
|и(г)-н А —>0, 
что и завершает доказательство сходимости линеаризованной задачи к решению исходной задачи. 
Теорема 2. Пусть дана начально-краевая задача: 





(1 г), 4» т затаив) в Цг=Рх(о,Т), Бек". 
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В 


Здесь К = } къ] < 0 сначальными условиями: 





Аб а 
оы Е зтааН 


НИИ... МИ и 
(ру-ро) 4) | ПФ 


Н(х, у,0) = Нь(х,у), Но(х,у)е С? (р)ос(р), ста у)Но Е с(Б), (ху)е р, 





условиями (4)-(8) на границе области р, Функция Н(х,у,!) — решение задачи класса "в (Цт)о с(Ц,) р 


та „)Н ЕС (Пт) : 
Пусть далее на временной сетке в. = а =и.т, = 0,1,....№М, М-т=Т } построена линеаризованная начально- 


краевая задача: 








(") 
ПыА = к("9. Тыс - зтааН(") а, , п=1,2,..., М. 
$1 фо 
Аба в № и 





Здесь #0) = С ‚ зтааН аа 1) 


- с начальными условиями: 
ут фо 








в "|1 
(р — Ро ) га) 
н(® (х, у.) = Но (х, у), н(®) (ху) — н("—) (ху), (х, у) Е Б, п=2..... М, 
граничными условиями (4)-(8), удовлетворяющими условию (17)-—(18). Функция НЯ (х, 5,1) — положительное реше- 


ние задачи класса С? ( НО < С ( Ц г) й &та4| ЭН ЕС (п ‚) . Также предположим, что существуют и непрерывны 
смешанные частные производные до третьего порядка включительно по пространственной и временной переменной # 


от функций Н (х, у, 1), Н(х, у,1) для (х,у)Е О, 0<Е<Т: 


д дн`\\) д ан`\\) д дн) д ОН 

— | этаа| — ||, —| этаа| — ||, —| етаа| — ||, —| этаа| — ||. 

ох 0 ду [6 дх 0 ду [011 

Тогда, если ("1 > К > 0, К =бомЕ К ЕС! (2), п=1,2,..., № , то имеет место оценка 


1 
НТ] < а” -[,-С-Т, (60) 


где С — постоянная величина, большая 0, т. е. положительное решение линеаризованной задачи Н ( х,у,Ё) стремится 
к решению Н (х, у,!) нелинейной задачи в норме Гильбертова пространства /, (2) при 0, М№-ч=Т. 


Заключение. Исследована линеаризованная двумерная математическая модель транспорта донных наносов в 
прибрежных системах, удовлетворяющая основным законам сохранения и учитывающая следующие физические про- 
цессы и параметры: пористость донного материала; критическое значение тангенциального напряжения, при котором 
начинается транспорт донных материалов; динамически изменяемая геометрия донной поверхности за счет движения 
водной среды. 

Получены условия положительности решений линеаризованной начально-краевой задачи и их сходимости к 
решению исходной нелинейной начально-краевой задачи в норме Гильбертова пространства Г, сформулированные в 
виде теорем. 
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